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ABSTRACT 

We report the discovery of a large, sudden, and persistent increase in the spin-down rate of B0540-69, 
a young pulsar in the Large Magellanic Cloud, using observations from the Swift and RXTE satellites. 
The relative increase in the spin-down rate i> of 36% is unprecedented for B0540-69. No accompanying 
change in the spin rate is seen, and no change is seen in the pulsed X-ray emission from B0540-69 
following the change in the spin-down rate. Such large relative changes in the spin-down rate are seen 
in the recently discovered class of “intermittent pulsars”, and we compare the properties of B0540-69 
to such pulsars. We consider possible changes in the magnetosphere of the pulsar that could cause 
such a large change in the spin-down rate. 

Subject headings: pulsars: individual (PSR B0540-69) 


1. INTRODUCTION 

B0540-69 is a young, rotationally powered pulsar in 
the Large Magellanic Cloud (LMC) that is similar to the 
Crab pulsar in many ways. With a rotation frequency 
1 / of almost 20 Hz, it is one of the most rapidly ro¬ 
tating young pulsars known. Its spin-down luminosity 
(—47r^Ji/f') of 1.5 X 10^®J45 erg s“^, where I 45 is the 
moment of inertia of the neutron star (NS) in units of 
10'^® g cm^, is also among the largest for all pulsars. The 
spin-down rate of a pulsar is often described in terms of 
a braking index n, in which = —kv^. The slowdown 
of B0540-69 is relatively stable for a young pulsar, and it 
is one of only 8 young pulsars for which a braking index 
has been reliably measured (Lyne et al. 2015). It has a 
characteristic spin-down age (—z//2z>) of ^ 1600 years. 

Pulsars are remarkably stable rotators, and the devi¬ 
ations from clocklike precision can provide information 
about the structure and processes at work in the neu¬ 
tron star and its magnetosphere. Sudden changes in ly 
and 1 / (“glitches”) are occasionally seen in pulsars, es¬ 
pecially those with ages between 10^ and 10® years (Yu 
et al. 2013). Glitches are thought to occur in rotationally 
powered pulsars when angular momentum is transferred 
from a more rapidly rotating component of the NS to the 
outer crust (e.g., Anderson & Itoh 1975; Eranco et al. 
2000). A different kind of rapid change is seen in some 
pulsars that affects the spin-down rate. In their exami¬ 
nation of long-term monitoring of 366 pulsars, Lyne et al. 
(2010) found that many pulsars rapidly switch between 
two different spin-down rates. The time scales for the 
transitions are quasi-periodic, with typical time scales of 
years. “Intermittent pulsars” are extreme examples of 
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such pulsars. They transition from a radio-on state to 
radio-off state with a simultaneous change in the spin- 
down rate. The best studied example is PSR B1931-I-24 
(Kramer et al. 2006; Young et al. 2013) for which mul¬ 
tiple transitions on a timescale of weeks from a radio 
bright state with i> = —16 x 10“^® Hz s“^ and a radio 
quiet state with = —10.8 x 10“^® Hz s“^ have been ob¬ 
served. Similar behavior on longer time scales has been 
reported for PSR J1832-I-0029 (Lorimer et al. 2012) and 
PSR J1841-0500 (Camilo et al. 2012). All three pul¬ 
sars show rapid transitions between two states with sta¬ 
ble spin-down rates, large differences in the spin-down 
rates, large changes in the radio flux, and no simultane¬ 
ous change in v. Unlike most explanations for glitches, 
models for these state changes have emphasized changes 
in the pulsar’s magnetosphere (e.g., Li et al. 2012). None 
of the intermittent pulsars have been detected in the op¬ 
tical, X-ray, or gamma-ray bands, but somewhat similar 
behavior has been seen for the gamma-ray pulsar PSR 
J2021-I-4026 (Allafort et al. 2013). A 4% increase in v 
was seen with a simultaneous decrease of 18% in the flux 
above 100 MeV with a time scale for the transition of 
less than a week. Since there is no radio counterpart, 
any change in radio flux is not known. 

B0540-69 has been extensively studied since its discov¬ 
ery with the Einstein Observatory (Seward et al. 1984). 
A small glitch for B0540-69 was reported by Zhang et al. 
(2001) with a relative change in 1 / of 1.9 x 10“® and h 
of 8.5 X 10“®. The reality of this glitch was disputed by 
Cusumano et al. (2003) and later supported by Living¬ 
stone et al. (2005). Ferdman et al. (2015) examined 15.8 
years of data from RXTE and report a second glitch with 
a relative change in i/ of 1.6 x 10“® and of 9.3 x 10“®. 
Both these changes in are orders of magnitude smaller 
than we report in this paper. No other glitches have been 
reported in the extensive monitoring of the pulsar. Op¬ 
tical pulsations are also seen (Mignani et al. 2010 and 
references therein). Manchester et al. (1993) discovered 
radio pulsations with an above average luminosity at 640 
MHz of 1200 mJy kpc^. 

In this paper we present and discuss new timing anal¬ 
ysis of B0540-69 as observed with RXTE and Swift. De¬ 
tails of the observations and the results of the temporal 
analysis are reported in Section 2, and interpretations 
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are given in Section 3. Finally, Section 4 is a summary. 

2. OBSERVATIONS AND DATA REDUCTION 

B0540-69 was observed with the Proportional Counter 
Array (PCA) on board the Rossi X-ray Timing Explorer 
(RXTE) (Bradt et al. 1993) during a 12.5-year campaign 
to monitor the nearby PSR J0537-6910. The final obser¬ 
vation was on Dec. 31, 2011. We concentrate on observa¬ 
tions covering the final 140 days of the campaign. Results 
from RXTE observations have been reported by Zhang 
et al. (2001), Livingstone et al. (2005) and Ferdman et al. 
(2015). The PCA is composed of five co-aligned xenon 
detectors (Proportional Counter Units (PCUs)) with a 
total area of ^ 6500 cm^. Individual PCUs were rou¬ 
tinely turned off and on to reduce the number of high- 
voltage breakdowns. Additional observations were made 
with the X-Ray Telescope (XRT) instrument (Burrows 
et al. 2005) on the Swift Gamma-Ray Burst Explorer 
(Gehrels et al. 2004) starting in Feb., 2015. Table [1] is a 
log of the observations. 

Data reduction for both missions followed very sim¬ 
ilar procedures. X-ray events were screened to maxi¬ 
mize the signal-to-noise ratio, and then photon arrival 
times were corrected to the solar system barycenter with 
the FTOOL® faxbary for RXTE and BARYCORR for 
Swift using the Hubble Space Telescope position of a = 
05/^40^11 202*, S = -69°19'54.17" (J2000.0) (Mignani 
et al. 2010). The mid-point of each observation was cho¬ 
sen as the epoch, and then the best period was deter¬ 
mined by comparing multiple folded light curves using 
EFSEARCH. Since the folded light curve is approximately 
a sine wave (Cusumano et al. 2003), the phase at the 
epoch was determined by fitting a sine wave to the folded 
light curve produced with efold and using the phase of 
the peak of the sine wave. The resulting frequencies and 
phases were fit to the usual truncated Taylor series ex¬ 
pansion of the phase and its time derivatives v, y, and 
i). 

Uncertainties quoted in this paper are given at the 
90% confidence level unless otherwise noted. For multi¬ 
parameter fits, the other parameters are allowed to vary 
when calculating the limits. 

2.1. RXTE/PCA Observations 

The PCA is sensitive to X-rays in the 2 to 60 keV band 
with moderate {AE/E ~ 18%) resolution. Each event 
is time-tagged on the spacecraft with an accuracy better 
than 5 x 10“^ s (Rots et al. 1998). To improve the signal- 
to-noise ratio, only events in the first xenon layer in the 
energy range of 3 — 20 keV were included. 

Figure [T] shows the frequency residuals for the last 
14 RXTE observations relative to the best-fit ephemeris 
model for the first 12 of the observations. The value of 
f is set to 3.249 x the best-fit value for obser¬ 

vations between Oct. 2010 and Dec. 2011. The model, 
whose parameters are given in Table [2J provides a very 
good fit to the first 12 observations with rms phase resid¬ 
uals of ^1%, but the final two observations deviate dra¬ 
matically. A linear fit to the frequencies for the final two 
observations indicates that has changed by—6.Ox 10“^^ 
Hz s“^ with a statistical uncertainty of 0.8 x 10“^^ Hz 
s“^ at some time between the observations on Dec. 3 and 

® http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 


Dec. 17. This is an increase of 32 ± 4% in the rate at 
which the pulsar’s spin is slowing. This value is consis¬ 
tent with, but less accurate than, the value determined 
in Section 2.2 using both Swift and RXTE data. We 
note that our value for is ^ 14% lower than the values 
reported by Ferdman et al. (2015) for longer intervals 
of RXTE data. Reprocessing the RXTE data using the 
Ferdman et al. values for i) demonstrates that revising i) 
would have a negligible effect on our results and would 
not change our conclusions. 

If the sudden change in v were due to a glitch, then 
a simultaneous change in v would be expected. Surveys 
of glitches in other pulsars find that the amplitude of a 
glitch in 1 / correlates with the amplitude of the glitch in 
v for pulsars in general (Espinoza et al. 2011) and for 
the Crab Pulsar in particular (Lyne et al. 2015). The 
correlation for pulsars in general would indicate an ac¬ 
companying glitch in B0540-69 larger than ~ 1 x 10“® 
Hz, and the correlation for the Crab would indicate an 
accompanying glitch of ^ 1 x 10“^ Hz. There is no in¬ 
dication of an accompanying change in z/, and we place 
limits on its possible size by assuming an instantaneous 
glitch with no recovery and extrapolating the fit to the 
post-glitch frequencies back toward the final pre-glitch 
observation. With these assumptions the largest glitch 
would occur immediately after the final pre-glitch obser¬ 
vation with a size of 1.2 x 10“® Hz. Later glitches would 
be smaller until at 200 ks after final pre-glitch observa¬ 
tion there would be no glitch in v. 

Large glitches in z> typically decay away on time scales 
from about a day to weeks. Since the new value for 
z> is consistent with that later seen with Swift (Section 
123), there is also no indication of such a recovery. We 
conclude that the sudden change in v is almost certainly a 
state transition rather than a glitch. We will refer to the 
state before the transition as the “low spin-down” state 
and the state after the transition as the “high spin-down” 
state. There is also no indication of a change in the 
pulsed X-ray emission from B0540-69. Figure [2] shows the 
average pulsed count rate {rpuise = A X C/texp/NpCU, 
where A is the relative sine wave amplitude, C is the 
total number of counts detected, t^xp is the exposure 
time, and Npcu is the average number of PCUs that 
were on during the observation). The average pulsed 
count rate for the two observations after the transition 
(0.733 ±0.086 s“^) is consistent with the average for the 
first 12 observations (0.718 ± 0.035 s“^). Further the 
folded light curves for these two observations show no 
significant change. 

To confirm the continued nominal performance of the 
PCA including the final two observations, we processed 
the same observations for PSR J0537-6910 following the 
procedures given in Marshall et al. (2004). After the 
glitch near MJD 55815, the final 12 observations are very 
well fit using the standard form for the ephemeris with 
an rms phase residual of ~ 1%. We conclude that the 
unexpected values of v for B0540-69 are not due to any 
instrumental effect. 

2.2. Swift XRT Observations 

The X-ray Telescope started observations of B0540-69 
on Feb. 17, 2015, as part of a Target of Opportunity 
campaign. The XRT is a focusing X-ray telescope with 
a CCD detector with an effective bandpass of 0.3 - 10 
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keV. All observations were made using the Window Tim¬ 
ing (WT) mode in which the central 200 CCD columns 
are continously read out, providing a time resolution of 
1.7 ms. The XRT data were reduced with the standard 
software (xrtpipeline vO.13.1) applying the default fil¬ 
tering and screening criteria (HEASoft 6.16), using 
the 20140709 update to the XRT CALDB files and the 
20150428 update to the clock correction file. The phase 
and frequency were measured for each continuous view¬ 
ing interval (a “snapshot”). Each of the Swift entries in 
Table [T] except the first consists of two snapshots sepa¬ 
rated by about 96 minutes (the duration of a Swift orbit). 
The first entry has a single snapshot. We obtained the 
best ephemeris by determining a model for the closely 
spaced observations and then refining the model as more 
widely separated observations were added. All the mod¬ 
els assumed a value of 3.249 x for based on 

measurements with RXTE before the state transition in 
2011. This parameter makes a very small contribution to 
the model because the duration of the current Swift cam¬ 
paign is only 65 days. The best-fit model was confirmed 
using a large grid search of possible combinations of p 
and b. The parameters of the best-fit model and the 90% 
confidence uncertainties are given in Tabled The model 
provides a good fit to the data with rms phase residuals 
of ~3%. Figure [U which shows the frequency residuals 
relative to the RXTE ephemeris before the transition, 
displays the overall history of ly. 

We used two methods to determine i> after the state 
transition. The first method uses the best-fit model for 
the XRT data (Table [5]) with a value of —2.52871 ± 
0.00008 X 10“^°s“^ at the XRT epoch. After adjust¬ 
ing by —3.33 X 10“^^s“^ to account for the effect of i>, 
this is a change in i/ from the value before the transition 
(Tabled!) of -0.66882 ± 0.00008 x lO-^^s-^. 

The second method compares the average value of p 
in the final two RXTE observations (19.72601295 Hz 
±4.6 X 10“® Hz at MJD 55919.71) with the best-fit value 
for the XRT data of 19.70077383 Hz ±4 x 10"® Hz. The 
resulting is —2.5230 ± 0.0005 x 10“^° s“^ at MJD 
56498.62. After adjusting by 1.62 x 10“^®s“^ to account 
for the effect of ii, the result is —2.5214 ± 0.0005 x 10“^° 
s“^ at the XRT epoch, which is larger than result from 
the first method by 7.3 ± 0.5 x 10“^® s“^. This discrep¬ 
ancy is shown in Figure [T] in which the extrapolation of 
the XRT ephemeris back to the final two RXTE obser¬ 
vations predicts a value for i/ that is 7.3 x 10“®s“^ larger 
than that observed. The uncertainty of the extrapolated 
values due to the uncertainty in the model’s v is only 
8 X 10“^ Hz. The discrepancy suggests that the evo¬ 
lution of the high spin-down state is more complicated 
than our simple model. The second method indicates 
a change in zz from the value before the transition of 
-0.6615 ±0.0005 X lO-^^s-^. 

Since the value of f after the state transition is uncer¬ 
tain and the statistical uncertainties are small, we use 
the value of the adjustments due to H to estimate the 
uncertainty. The relative change in b is then 35.9 ±0.2% 
and 35.5 ± 0.1% for the first and second methods respec¬ 
tively. We adopt 35.7 ± 0.4%, the average of the two 
methods with an uncertainty that encompasses the in¬ 
dividual error bars, as our best estimate of the relative 
change in z>. 

The folded light curve is similar to that reported from 


other X-ray observations (e.g., Cusumano et al. 2003). 
Although the actual light curve is more complicated, a 
model of a constant plus a Gaussian provides a good de¬ 
scription of the XRT data. The best-fit standard devia¬ 
tion for the Gaussian is 0.174±0.022, which is consistent 
with a fit to the folded light curve for the RXTE data 
before the state transition. 

3. DISCUSSION 

B0540-69 has been extensively monitored by numer¬ 
ous observatories in the more than 30 years since its dis¬ 
covery including Ginga observations spanning 4.4 years 
(Deeter et al. 1999) and the 12.5 years of observations 
with RXTE. More than 10 measurements of the braking 
index for the pulsar, which require long spans of data, 
have been reported (Cusumano et al. 2003, who list pre¬ 
vious measurements; Ferdman et al. 2015). Values for 
the braking index range from 1.81 (Zhang et al. 2001) to 
2.74 (Ogelman & Hasinger 1990). 

None of these studies found a sudden change in z> com¬ 
parable to the one reported here. The glitches in v re¬ 
ported by Zhang et al. (2001) and Ferdman et al. (2015) 
are ~2000 times smaller, and knowledge of the pulsar’s 
phase is maintained throughout the 12.5 years of RXTE 
monitoring until the final two observations. When dis¬ 
covered with the Einstein Observatory (Seward et al. 
1984), B0540-69 had a z> of -1.900 x 10-i° Hz s-\ which 
indicates that the pulsar was in the low spin-down state 
in 1979 and 1980. The lack of published values of f' in¬ 
dicative of the high spin-down state suggests that B0540- 
69 remained in the low spin-down state the vast majority 
of the time, if not all the time, until the transition in late 
2011. 

The sudden change in the spin-down rate is most likely 
associated with a change in the magnetosphere of the 
NS. One possibility is a global change in the conduc¬ 
tivity of the magnetosphere. Li et al. (2012) developed 
solutions for pulsar magnetospheres with finite resistiv¬ 
ity that mitigate some of the limitations of the vacuum 
and force-free models. The dependence of the spin-down 
luminosity was calculated for a wide range of the con¬ 
ductivity parameter a (expressed in terms of the pulsar 
angular rotation rate D) and inclination angle a using 
a three-dimensional numerical code. For broad ranges 
of a (0° to 90°) and the conductivity ((ct/H)^ from 0.04 
to 4.0), a 36% increase in the luminosity (and by exten¬ 
sion in I z/ I) can be achieved with an increase in a of less 
than a factor of 6. There are maximum values of (cr/D)^, 
ranging from 40 at a = 0° to 1.3 at a = 90° above which 
it is no longer possible to explain an increase in the lu¬ 
minosity with an increase in cr. Since for large values of 
(cr/D), the component of the electric field parallel to the 
magnetic field E\\ oc (Kalapotharakos et al. 2014), 
an increase in cr may reduce the gamma-ray flux from 
the pulsar. The long-term light curve of B0540-69 in the 
energy range of 200 MeV to 100 GeV (Fermi-LAT Col¬ 
laboration 2015) shows no indication of change in the 
gamma-ray flux with an upper limit of 30%. Detailed 
modelling, such as that done by Brambilla et al. (2015) 
for other Fermi pulsars, is needed to understand the im¬ 
plication of this constraint on parameters for B0540-69, 
but this is beyond the scope of this paper. 

Another possible explanation for the spin-down transi¬ 
tion is a change in torque due to a change in the plasma 
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outflow from the polar caps. Kramer et al. (2006) sug¬ 
gested this to explain the state transitions in the inter¬ 
mittent pulsar PSR B1931-f24. Unlike PSR B1931-k24, 
B0540-69 appears to be radio bright in its low spin-down 
state, but its radio flux may significantly increase or its 
pulse shape may change in the high spin-down state. Fol¬ 
lowing Young et al. we estimate the change in the charge 
density Ppiasma as 7.1 x C m“^, assum¬ 

ing a pulsar radius of 10 km and I 45 of 1. This estimate of 
0.78 C m“^ is close to the Goldreich-Julian charge den¬ 
sity for B0540-69 of 1.09 C m“^. An increase in Ppiasma 
predicts an increase in the radio flux, which can be tested 
with radio observations of the high spin-down state. It 
also predicts a reduction in the braking index as the rel¬ 
ative importance of magnetic dipole radiation decreases. 

A change in the amount of open poloidal magnetic flux 
would also change the spin-down luminosity, which is ex¬ 
pected to be proportional to where Topen is the ra¬ 
dius beyond which the magnetic field lines become open 
(Contopoulos 2007). An 36% increase in the spin-down 
luminosity with no change in v would require a 15% de¬ 
crease in Topen- 

4. SUMMARY 
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The spin-down rate of B0540-69 increased suddenly by 
36% (a change in z> of —6.7 x 10“^^ Hz s“^) between Dec. 
3 and Dec. 17, 2011. Observations with Swift in 2015 
show that this change in spin-down rate has persisted 
for more than 3 years. Such a large persistent change in 
the spin-down rate has never been reported for B0540- 
69. There is no indication of a simultaneous change in the 
pulse shape in X-rays, the pulsed X-ray luminosity, or the 
total gamma-ray luminosity. We interpret this change 
as a transition between two stable states of the pulsar, 
similar to state transitions in other pulsars. B0540-69 
extends the class of state changing pulsars to include a 
very young and luminous example. The transitions in 
B0540-69 appear to be rare with the low spin-down state 
probably lasting more than 30 years. Future work can 
test predicted changes in the radio flux and the braking 
index in the high spin-down state. 
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Table 1. Observing Log 


Satellite 

Ohs. ID 

Date 

(1) 

MJD 

(1) 

Texp 

(2) 

RXTE 

96023-01-17 

2011 Aug 13 

55786 

6.6 

RXTE 

96023-01-42 

2011 Sep 06 

55810 

7.3 

RXTE 

96023-01-43 

2011 Sep 14 

55818 

7.5 

RXTE 

96023-01-44 

2011 Sep 18 

55822 

6.9 

RXTE 

96023-01-20 

2011 Sep 23 

55827 

6.9 

RXTE 

96023-01-45 

2011 Sep 26 

55830 

7.2 

RXTE 

96023-01-46 

2011 Sep 26 

55830 

6.3 

RXTE 

96023-01-21 

2011 Oct 09 

55843 

6.7 

RXTE 

96023-01-22 

2011 Oct 22 

55856 

6.9 

RXTE 

96023-01-23 

2011 Nov 04 

55869 

6.9 

RXTE 

96023-01-24 

2011 Nov 18 

55883 

7.3 

RXTE 

96023-01-25 

2011 Dec 03 

55898 

6.6 

RXTE 

96023-01-26 

2011 Dec 17 

55912 

7.1 

RXTE 

96023-01-19 

2011 Dec 31 

55926 

7.8 

Swift 

00033603002 

2015 Feb 17 

57070 

1.2 

Swift 

00033603004 

2015 Feb 23 

57077 

2.1 

Swift 

00033603005 

2015 Feb 25 

57078 

1.8 

Swift 

00033603006 

2015 Feb 25 

57078 

2.4 

Swift 

00033603007 

2015 Mar 11 

57092 

2.1 

Swift 

00033603008 

2015 Apr 11 

57123 

2.3 

Swift 

00033603009 

2015 Apr 13 

57125 

2.3 

Swift 

00033603010 

2015 Apr 23 

57135 

2.5 


Note. — (1) At the start of the observation; (2) Exposure time in ks 


Table 2. Ephemeris Parameters 


Parameter 

Pre-transition 2011 
(1) 

2015 

(1) 

Epoch (Modified Julian Date) 

55892.352587464 

57077.537772223 

Phase. 

0.000 (15) 

0.000 (23) 

(Hz). 

19.72655182 (1) 

19.70077383 (4) 

V (10“l°s-2). 

-1.86322 (2) 

-2.52871 (8) 

V (10-2ls-3). 

3.249 (fixed) 

3.249 (fixed) 


Note. — (1) The numbers in parentheses are the 90% confidence errors quoted in the last digit. 
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Fig. 1.— The frequency residuals for the RXTE and XRT observations relative to the ephemeris model (the dark blue dashed line) 
determined with the RXTE observations before the state transition. The dash-dotted red line shows the best-fit ephemeris for the XRT 
data. One-cr uncertainties are shown, but they are smaller than the symbols for the RXTE data. 
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Fig. 2.— The pulsed count rate with 1-cr uncertainties for the RXTE observations. 






















